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Abstract. The diffuse continuum emission from the Galactic plane in the energy range 18-1000 keV has been 
studied using 16 Ms of data from the SPI instrument on INTEGRAL. With such an exposure we can exploit 
the imaging properties of SPI to achieve a good separation of point sources from the various diffuse components. 
Using a candidate-source catalogue derived with IBIS on INTEGRAL and a number of sky distribution models 
we obtained spectra resolved in Galactic longitude. We can identify spectral components of a diffuse continuum of 
power law shape with index about 1.7, a positron annihilation component with a continuum from positronium and 
the line at 511 keV, and a second, roughly power-law component from detected point sources. Our analysis confirms 
the concentration of positron annihilation emission in the inner region < 10°), the disk (10° < \l\ < 30°) being 
at least a factor 7 weaker in this emission. The power-law component in contrast drops by only a factor 2, showing 
a quite different longitude distribution and spatial origin. Detectable sources constitute about 90% of the total 
Galactic emission between 20 and 60 keV, but have a steeper spectrum than the diffuse emission, their contribution 
to the total emission dropping rapidly to a small fraction at higher energies. The spectrum of diffuse emission is 
compatible with RXTE and COMPTEL at lower and higher energies respectively. In the SPI energy range the 
flux is lower than found by OSSE, probably due to the more complete accounting for sources by SPI. The power- 
law emission is difficult to explain as of interstellar origin, inverse Compton giving at most 10%, and instead a 
population of unresolved point sources is proposed as a possible origin, AXPs with their spectra hardening above 
100 keV being plausible candidates. We present a broadband spectrum of the Galactic emission from 10 keV to 
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100 GeV. 
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1. Introduction 




Interstellar gamma-ray emission from the Galaxy is a ma- 
jor study objective for INTEGRAL. The inner Galactic 
ridge is known to be an intense source of continuum hard 
X- and soft 7-r ay emission: hard X -ray emission was dis- 
covered in 1972 l|Bleach et all 1 197^) . and interstellar emis- 
sion has subsequently been observed from keV to MeV en- 
ergies by ASCA, Ginga, RXTE, OSSE, COMPTEL and 
most recently by Chandra and XMM-Newton. Most di- 
rectly comparable to our INTEGRAL analysis are results 
from OSSE on the Compton Gamma Ray O bservatory 




iDogiel et al At MeV energies the o rigin of the 



Continuum emission of diffuse, interstellar nature is 
expected in the hard-X and 7-ray regime from the physi- 
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Alternatively, the origin of the ridge emission could 
be attributed to a population of sources too weak to 
be detected individually, and hence it could be not 
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truly interstellar. Gamma-ray telescopes in general have 
inadequate spatial resolution to clarify this issue. For 
X-rays, important progress in this are a was made by 
ASCA (jKaneda et all Il997j) and Ginga l)Yamasaki et all 
119971) . More recently, high-r esolution imag i ng in X-rays 
(2-10 keV) with Chandra l|Ebisawa et all 1200 ll 12005^ 
has proven the existence of a truly diffuse compo- 
nent; similarly, it h as been inferred from XMM-Newton 
l|Hands et all l2004j) that 80% of the Galactic-ridge X- 
ray emission is probably diffuse, and only 9% can be ac- 
counted for by Galactic sources (the rest being extragalac- 
tic sources). 

The study of the Galactic-ridge continuum emission is 
a key goal of the INTEGRAL mission. The high spectral 
resolution combined with its imaging capabilities promises 
new insights into the nature of this enigmatic radiation. 
Furthermore, reliable modelling of the diffuse emission will 
be essential for the study of point sources in the inner 
Galaxy, since it constitutes a large and anisotropic back- 
ground against which sources must be resolved. 

Previous work based on first, smaller sets of 
INTEGRAL/SPI observations have reported the detec- 
tion of diffus e emission at a level consistent wit h previous 
experiments ^Strong et all l2003t IStrond. l2003^> . But sta- 
tistical and systematic errors were large, due in part to the 
uncertainty in the point-source contributio n. Meanwhile , 
a new analysis of INTEG RAL /IBIS data l|Lebrun et all 
12004 iTerrieret alllioO^ showed that, up to 100 keV, in- 
deed a large fraction of th e total emissi o n from the inner 
Galaxy is due to sources. IStrong et alJ (|2004ch used the 
source catalogue from this work (containing 91 sources) 
as input to SPI model fitting, giving a much more solid 
basis for the contribution of point sources in such an analy- 
sis. This exploited the complementarity of the instruments 
on INTEGRAL for the first time in the context of diffuse 
emission. In the present work we extend this analysis to 
include a much larger observation database, and use the 
2nd IBIS Catalogue to provide source candidates. 

We also d raw attention to a related paper 
ijBouchet et all 120051) . which describes an indepen- 
dent but similar study of SPI data, but focussing more 
on the point-source contribution. 



We use data from the SPI (INTEGRAL Spectrometer) 
instrument; descripti ons of the inst r ument and perfor- 
mance are given m IVedrenne etaH l|2003|) : lAttie et alJ 
(2003). The energy range covered by SPI is 20 keV - 
8 MeV, but here we restrict the analysis to energies up 
to 1 MeV; above these energies the statistics are small 
and the analysis is more difficult, so is reserved for future 
work. The data were pre-processed using the INTEGRAL 
Science Data Centre (ISDC) Standard Analysis software 
(OSA) up to the level containing binned events, with 0.5 
keV bins, with pointing and livetime information. GCDE 
data from about 200 orbital revolutions from 15 - 259 
were used. 8100 pointings were used; the sky exposure is 
shown in Fig^ the total exposure livetime is 1.6 x 10 7 s. 
The exposure per pointing is typically 1800 s. The energy 
calibration is performed using instrumental background 
lines with known energies; while this is a critical opera- 
tion for line studies (where sub-keV accuracy is required), 
for continuum studies a standard calibration (~1 keV ac- 
curacy) is quite adequate. Various energy binnings were 
used, depending on the available statistics as a function 
of energy. Only single-detector events are used here, since 
they dominate below 1 MeV. After the failure of detector 
2 at LID 1 1434 the remaining 18 detectors were used, and 
after the failure of detector 17 at IJD 1659 the remaining 
17 were used. 

The instrumental response is based on e xtensive M onte 
Carlo simulations and parameterization ijSturner et all 
2003); this has been tested on the Crab in-flight calibra- 
tion observations and shown to be reliable to better than 
20% in absolute flux at the current state of the analysis. 

Since INTEGRAL data are dominated by instrumental 
background, the analysis needs to have good background 
treatment methods. In the present work the background 
ratios between detectors are obtained by averaging the 
entire dataset over time; this has been found to give results 
which hardly differ from taking 'OFF' observations, and 
has the advantage of much higher statistics. A check on 
the systemat i cs in t his background approach was given in 
IStrong et al\ l|2003j) . 



2. INTEGRAL Observations 

The INTEGRAL Core Program ()Winkler et all 12003^ for 
A01 includes the Galactic Centre Deep Exposure (GCDE) 
which maps the inner Galaxy (330° < I < 30°, -20° < b < 
20°) with a viewing time of about 4Ms per year. The full 
region is covered in one GCDE cycle, and there are two 
cycles per year. Subsequent AOs have varied this strategy 
somewhat. In addition the Galactic Plane Survey (GPS) 
covers the whole plane at low Galactic latitude, with lower 
exposure than the GCDE. Data from the first three years, 
including GCDE, GPS, data which has become public up 
to March 2005, and Open Time data for which permis- 
sion has been obtained from the corresponding Principal 
Investigator, are used for the study reported here. 



3. Overview of Methods 

3.1. Data combination, energy binning, selection 

The program spiselectscw combines the observation-by- 
observation standard-processed data into a single set of 
datasets suitable for input to analysis software. Its func- 
tions include rebinning in energy, rejection of anomalous 
data (e.g. high rates due to solar flares), and detector se- 
lection (singles, multiples). This step is required to reduce 
the enormous data volume of the standard-processed data 
(with its 0.5 keV binning) into a manageable form. 



IJD=INTEGRAL Julian Date starts at 1 Jan 2000 
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3.2. Background template generation 

The program spioffback generates a background template 
by averaging the detector ratios over all pointings (either 
from a set of high-latitude 'OFF' observations or from the 
ON observations themselves), and using these ratios to 
compute predicted counts using the livetime information. 
The total predicted counts are normalized to the total ob- 
served counts in each energy interval: this normalization 
is for convenience, but is not essential since the absolute 
background level is determined later in the fitting proce- 
dure. The template is generated separately for each of the 
periods separated by detector 2 and 17 failures, since these 
failures have a signficant impact on the detector ratios. 

3.3. Model fitting 

The program spimodfit uses the binned count data, point- 
ing and exposure information, and the instrumental re- 
sponse, to fit sky models, point sources and background. 
The output is in the form of FITS files containing, for all 
processed energies: the fit parameters, the covariance ma- 
trices and the input models, spimodfit is a development 
(by H. Hallo in, MPE) from spidiffit w hich was used in 
earlier work I Strong et all 120031 l2f)04ch . It includes more 
flexible background fitting, optimized fitting routines and 
memory use, and improvement formatting of the output. 
We use spimodfit version 2.4 in this work. 

The model skymaps to be fitted are generated with the 
program gensky, which includes HI and CO survey data, 
as well as Gaussians of arbitrary position and width. The 
input source catalogue is in the ISDC catalogue format. 



mainly on the level of the power-law component; for the 
positronium the effect is small. 

3.5. Evaluation and plotting 

The program galplot combines the results from many spi- 
modfit runs into spectra integrated over sky regions, and 
sources. It evaluates the fluxes and error bars using the fit 
parameters covariance matrices. It plots the SPI spectra 
together with results from other instruments like RXTE, 
COMPTEL, E GRET, and t heoretica l models from the 
galprop project IStrong et alJ lj20f)4; 



4. Model Fitting 

4.1. Basic approach 

We distinguish image space and data space in the usual 
way, and define the instrument response as the relation 
between them. The image is Ij and the expected data is 
dk- The expected background is bk Let Rjk be the response 
of data element k to image element j. Then 



dk — 2_, Rjkl] + bk 



For the the special case of an image component consist- 
ing of a point source, only the position is required and 
computing dk requires no convolution. 

An image model is a parameterized algorithm for com- 
posing an image from components. For a linear model, 



3.4. Spectra 

The appropriate size of the energy bins is a trade-off 
between the high instrumental resolution and the avail- 
able statistics. Since this cannot be predicted a priori, 
we make the fits in a range of binsizes over the whole 
spectrum, and choose for presentation those for which a 
satisfactory signal-to-noise is achieved. Hence the binsize 
increases from ~10 keV at low energies to ~100 keV at 
high energies. 

The spectral fitting is done using the full energy re- 
distribution matrix in the form of 'RMF' and 'IRF' SPI 
response files (see Section 2 for references). The proce- 
dure used is rigorous for the case where the ratio between 
the three components of the IRF (photopeak and the two 
Compton components) is independent of incidence direc- 
tion: this variation is small 2 so the procedure is accu- 
rate. The energy response has a significant effect on the 
fits as can be seen from the difference between the raw 
data points and the fitted spectra (FigsEHSJ- The effect is 



2 The rms deviation of these component ratios with direc- 
tion, weighted by the sensitive area in that direction, is at 
the few percent level relative to the total. For example at 500 
keV: photopeak: 0.675 ± 0.025, Compton-1: 0.124 ± 0.0066 , 
Compton-2: 0.202 ± 0.03. 



where 9i are the model parameters. 

More generally, the image will still be described by a 
sum of components, but the image components will be 
non-linear functions of the parameters (e.g. Gaussian, ex- 
ponential) and each component is described by several pa- 
rameters 6i\ 



E 



Similarly the background can be constructed from 
components of a background model Bik 



bk = £ OiBik 



where 9i now introduces background parameters. The 
sums in the above expressions are over the appropriate 
subsets of parameters for image and background model re- 
spectively. In this way we can treat image and background 
model in the same way in the subsequent analysis, and 9i 
includes both. The only formal difference between image 
and background model is that the image is convolved with 
Rjk and the background is not: 



i=Ni 
j i=l 



1-Ni+Nb 
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Table 1. Diffuse model components used in the fitting 



Component 


description 


longitude 


latitude 






FWHM (deg) 


FWHM (deg) 


1 


HI 






2 


CO 






3 


Gaussian 


40 


5 


4 


Gaussian 


60 


5 


5 


Gaussian 


80 


5 


6 


Gaussian 


5 


5 


7 


Gaussian 


10 


5 


8 


Gaussian 


20 


5 


9 


Gaussian 


10 


10 



where there are Nj image components and Nb back- 
ground components, and N = Nj + Nb- 

In our modelling approach, the measured signal is rep- 
resented through model components: 

dk = ^2 diSik 

i 

where the sky part (diffuse + sources) is: 

S ik =J2 R okMij, i = l,Nj 

3 

and the instrumental-background part is: 

S lk = B lk , i = Nj + l,Nj + N B 

Details of the statistical model fitting method are given 
in Appendix A. 

5. The models 

5.1. Diffuse Galactic emission 

Model maps to represent the emission from cosmic-ray 
interactions wi th the gas a r e integ rated HI and 12 CO 
as described in IStrong et alJ l)2004a|i . Other components, 
including positronium, positron line, bulge and Galactic 
plane, are represented by a set of Gaussians of various 
widths summarized in Table 1. For the positronium and 
posi tron line we include Gaussians wit h FHWM = 5°, 
10° l|Jean et all 120031: iKnodlseder et alll2003[) . Generally 
only the sum of the components are presented since the 
separate components do not necessarily have physical 
significance; the idea is to include sufficient aspects of 
Galactic structure to be able to reproduce the total emis- 
sion with reasonable confidence. Summing subsets of com- 
ponents allows a morphological decomposition of the spec- 
trum for example into disc and bulge emission. 

Note that our model components are not orthogonal, 
but this is not a problem since our use of the covariance 
matrix explicitly handles the parameter correlations and 
the computation of the errors on linear combinations of 
components. 



5.2. Sources 

The sources fitted are b ased on a pr eliminary version of the 
2nd ISGRI Catalogue ijBird et all 1200.^ . containing 209 
sources. 3 Only the source coordinates are used here, since 
the flux in each energy range is determined in the fitting 
procedure simultaneously with the other components. 

At low energies it is known that most of the catalogue 
sources are detectable and make a corresponding contri- 
bution to the emission. At high energies on the other 
hand, only a few sources have a hard enough spectrum 
to be of relevance. Including the whole catalogue in fact 
leads to a 'glow' or 'cross-talk' phenomenon, whereby the 
diffuse emission is attributed to the sum of many cat- 
alogue sources at a very low level (< la), and the fit 
cannot distinguish between this source glow and emission 
from diffuse components. As a result in this case we get 
large error bars on both components, with highly anti- 
correlated fluxes, while the total flux (diffuse + sources) 
is still well constrained. Hence at higher energies the num- 
ber of sources included must be reduced using our knowl- 
edge of the individual sources from IBIS, SPI and other 
missions. Beyond 200 keV only a few sources are required. 
A systematic treatment of this effect is to include only 
sources above a given detection threshold in each energy 
range. We use a 3a threshold, thereby including all sources 
which could contribute while excluding the low-level 'glow' 
sources. Reducing this threshold (to la, 2a) allows an es- 
timate of the sensitivity of the diffuse spectrum to the 
chosen level. The energy bands below 100 keV are only 10 
keV wide (to exploit the energy resolution of SPI as far as 
allowed by statistics) , but use of these bands for the source 
selection is not sensitive enough and leads to extra fluctu- 
ations of the number of sources with energy; hence we use 
wider bands of 30 keV for this purpose. This reduces such 
fluctuations to an acceptable level.. Between 100 and 200 
keV the source-selection bands are widened from 40 keV 
to 60 keV, above 200 keV no widening is required. 

Table 2 shows the number of sources above a range of 
<7-levels, as a function of energy. 

Galactic sources are generally variable, except for some 
like the Crab nebula. At high energies where the sources 
are anyway unimportant, there is no problem, but at low 
energies the variability must be accounted for. The present 
analysis can take this into account by allowing the fluxes 
to vary on some timescale, with a corresponding increase 
in the number of free parameters. Only for strong sources 
is the variability expected to have an influence on the 
diffuse emission parameters. We investigate this by com- 



An analysis was also made with version 20 of the ISDC 
reference catalogue (available from http://isdc.unige.ch), in 
which sources detected by IBIS, SPI and JEMX are flagged; 
This catalogue includes published sources and those announced 
ATEL/IAU telegrams. There are 203 sources flagged as de- 
tected by IBIS in this catalogue. The results for the diffuse 
emission were not significantly different to those using the 2nd 
ISGRI Catalogue. 
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Fig. 1. Sky distribution of exposure, at 1 MeV. Galactic 
coordinates centred on I = b = 0. Units: 10 7 cm 2 s. Green: 
0.5-1.4, yellow: 1.4-4, orange: 4-13, red: 13-35 

paring the results assuming constant sources and various 
variability timescales. 

Table B.l shows the effect of the source variability on 
the determination of diffuse emission, for various energies 
and variability timescales. An important conclusion is that 
the effect is small, in fact it is within the quoted statis- 
tical errors for the case of a constant flux; so assuming 
constant source fluxes has no significant influence on the 
fitted diffuse intensity. Even al lowing for vari a bility of a 
few strong sources (taken from lBouchet et all J2005)) on 
the SPI pointing timcscale only changes the diffuse emis- 
sion results at the 20% level, and that only at the low- 
est energies. This is not really surprising since the diffuse 
emission has a very different signature in the data from 
any given source, and while an inaccurate source flux at a 
given epoch degrades the fit for that period, it has hardly 
any effect on the global fit. Allowing for variability of all 
sources on timescales above 50 days tends to reduce the 
diffuse flux slightly at low energies, by 50% in the most 
extreme case in the lowest energy range, otherwise by less 
than 30%, and always within the quoted errors. Above 
40 keV there is no effect. 

6. Results 

6.1. Background 

Fig [3 shows the relative background variation determined 
in the fitting, per pointing, for a typical energy range; the 
variations are as expected from tracers such as saturated 
Germanium rates and the general upward trend due to the 
solar cycle. Discontinuities appear at the times of detector 
2 and 17 failures, as expected. Occasional high values are 
encountered which were not high enough to be removed 
by the preliminary data filtering, but these are accounted 
for in the fitting procedure. 

6.2. Spectra 

Fig-Glshows the spectra of the summed catalogue sources 
and of total diffuse emission, for the inner region of the 




A/ 



! if 



Fig. 2. Time dependence of background for 48 - 58 keV de- 
termined from the fitting procedure. Relative units. Time 
axis in days since 1 Jan 2000 
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Fig. 3. Spectrum of diffuse emission summed over com- 
ponents (blue), and sum of all sources (red), for the inner 
region of the Galaxy (330° < I < 30°, |6| < 10°). 3cr source 
threshold. 



Galaxy (330° < I < 30°, |b| < 10°). We find a clear change 
from source-dominated to diffuse-dominated emission in 
the inner Galaxy, going from low to high energies. The 
transition energy is in the vicinity of 100 keV, but this 
depends on the fitted sky region due to variation of rela- 
tive intensities of the components along the plane of the 
Galaxy. The 'diffuse' emission may of course include con- 
tributions from unresolved sources. 4 

We have studied the effect of reducing the threshold 
for sources included in the fit from 3a to 2a and la. As 
expected the diffuse component decreases while the source 
component increases. However the diffuse component de- 



4 In this paper we refer to unresolved emission generically as 
'diffuse', independent of its true nature. 



6 Strong A.W. et al.: Gamma-ray continuum emission from the inner Galaxy 

Table 2. Number of sources detected above various significance levels, as function of energy, for the entire input 
catalogue. 



Energy 
(keV) 


1(7 


2cr 


3(7 


4(7 


5(7 


10ct 


20(7 


50a 


100(7 


18 0- 48 


145 


134 


123 


114 


93 


59 


33 


15 


8 


28 0- 58 


145 


123 


97 


82 


71 


39 


20 


g 


(j 


38 0- 68 n 


123 


84 


62 


47 


40 


17 


12 


4 


2 


48 0- 78 f) 


108 


63 


44 


34 


26 


14 


6 


3 


2 


58 0- 88 n 


100 


52 


36 


29 


22 


11 


6 


2 


2 


68.0- 98.0 


84 


57 


41 


27 


19 


12 


7 


2 


2 


78.0- 108.0 


77 


48 


30 


20 


18 


9 


5 


2 


2 


88.0- 148.0 


73 


38 


23 


17 


14 


8 


3 


2 


1 


128.0- 188.0 


48 


20 


12 


8 


6 


3 


2 


1 


1 


178.0- 258.0 


41 


12 


4 


3 


3 


2 


2 


1 





258.0- 338.0 


43 


9 


4 


2 


2 


2 


1 


1 





338.0- 418.0 


31 


4 


3 


2 


1 


1 


1 








418.0- 498.0 


25 


3 


1 


1 


1 


1 













Fig. 4. Spectrum of diffuse emission summed over model 
components, for the inner region of the Galaxy (330° < 
/ < 30°, |6| < 10°). 3cr source threshold. Also shown is 
the fit of SPI data to power-law plus positronium plus 
line, using the full energy response. Note that the fit lies 
below the data points because of the non-diagonal energy 
reponse of SPI. Magenta point: diffu se emission measured 
by IBIS/ISGRI llTerrier et al.ll2004h . 



creases by only 30% for the 2a source selection, which is at 
the lower limit for a plausible cut. For the lcr source selec- 
tion, the source component shows contamination from the 
diffuse emission above 100 keV due to 'cross-talk' between 
the components. We adopt the 3a selection as standard in 
what follows. 

The spectra can be fitted by a combination of power- 
law and the positron annihilation components (positro- 
nium continuum with the sharp upper edge at 511 keV 
and a linear decrease towards lower energies, and the 511 



Fig. 5. Spectrum of diffuse emission summed over model 
components, for the central part of the Galaxy (350° < 
I < 10°, |6| < 10°). See Fig 4 caption for details of fit. 



keV line). Fit results for the inner region of the Galaxy 
(330° < I < 30°, |6| < 10°), the central part of the Galaxy 
(350° < I < 10°, |&| < 10°), and the Galactic-disk re- 
gion (330° < I < 350°, 10° < I < 30°, |6| < 10°) are 
given in Table 3. Fig. 01 shows the spectrum in the inner 
region of the Galaxy (330° < I < 30°, |6| < 10°), and 
Fig. the spectrum for the central part of the Galaxy 
(350° < I < 10°, \b\ < 10°). All regions show a hard power 
law, with spectral index around 1.7. 

The positronium edge at 511 keV is highly significant, 
and below it the positronium continuum with its char- 
acteristic linear increase (I=AE) is clearly visible. For the 
Galactic-disk region (330° < I < 350°, 10° < I < 30°, \b\ < 
10°), Fig. the positronium component is small and not 
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Table 3. SPI diffuse emission intensity spectrum fitted to 
1(E) = A E"T + C E (E<0.511) + D<S(E-0.511) cm" 2 
sr" " 1 s _1 McV -1 where E is energy in MeV. 



energy, MeV 



Fig. 6. Spectrum of diffuse emission summed over model 
components, for the Galactic-disk region (330° < I < 
350°, 10° < I < 30°, |&| < 10°). See Fig 4 caption for 
details of fit. 



significant, as expected, and the spectrum is consistent 
with a power law. 

Using the fit results for the central part of the Galaxy 
(350° < I < 10°,\b\ < 10°), the 511 keV flux is 
0.81 ± 0.07 x 10 -3 cm~ 2 s _1 , and the 37 continuum flux is 
4.15±0.64 x 10~ 3 cm _2 s _1 . This gives a positronium frac- 
tion 5 / of 0.97—1.08 , consistent w i th oth e r determinations 
llKinzer et all l200lt Ijean et all 120041: IChurazov et~afl 
l2005lk although this is not the main goal of this analy- 
sis, it provides a consistency check on the analysis meth- 
ods. While the positronium flux decreases by a factor at 
least ^7 from \l\ < 10° to |/| > 10°, the power-law com- 
ponent decreases by only a factor ~ 2. This clearly shows 
the difference in angular distribution and hence a different 
spatial origin of the emission. 

The power-law continuum flux in this inner region is 
1.3 ± 0.25 x 10~ 3 cm- 2 s~ 1 for 300-500 keV, i.e. about 
1/3 of the positronium continuum in this range. The 
positronium component is small in the Galactic-disk re- 
gion (330° < I < 350°, 10° < I < 30°, |6| < 10°), which 
allows a more precise determination of the non-positron 
continuum. 

For comparison t he low-energy diffus e emission detec- 
tion by IBIS/ISGRI l)Terrier et all 12004^ is also shown in 
Fig El; the agreement is satisfactory. The SPI spectrum is 
compared with RXTE, COMPTEL and OSSE results in 
Fig H ; the comparison is made in the region the central 
part of the Galaxy (350° < Z < 10°, |fo| < 10°), since this 
is most compatibl e with the r egions reported for RXTE 
llRevnivtsevl. l2003h and OSSE llKinzer et allll999l) . 



parameter 




value 




error 




330° < I < 30°, 


b\ < 10° 










A 




8.32 10" 




1.8 10" 




7 




1.69 




0.09 




C 




1.19 10" 


1 


2.8 10" 


-2 


D 




2.36 10" 


-3 


3.8 10" 


-4 


350° < I < 10°, 


b\ < 10° 










A 




1.15 10" 


T~ 


2.2 10" 


-3 


7 




1.82 




0.08 




c 




2.61 10" 


1 


4.0 10" 


-2 


D 




6.62 10" 


3 


5.6 10" 


-4 


330° < / < 350° 


10° < I < 30° 










\b\ < 10° 












A 




7.87 L0 


~r- 


2.5 10" 




7 




1.47 




0.14 




C 




4.49 10" 


2 


3.2 10" 


-2 


D 




2.70 10" 


4 


4.3 10" 


-4 



7. Discussion 

7.1. Spectral shape 

The OSSE spectrum is somewhat higher than the SPI 
result in the 20-100 keV range (FigUJ. This suggests that 
the OSSE result includes significant source contamination, 
which is better accounted for in the SPI analysis. 

The RXTE spectrum is certainly above the extrapola- 
tion of the SPI power-law. This suggests a much steeper 
spectrum below 20 keV, or an additional component with 
a cutoff around this energy. For RXTE, we use the results 
of iRevnivtsel l|2003h 6 since this is pr obably more free of 
point- source contamination than e.g. IValinia fc Marshall 
(1998) and IValinia et al. I ll2000alk however it refers to 
the region 350° < I < 10°, 2 < \b\ < 10° , i.e. avoid- 
ing the Galactic ridge itself. The latitude distribution 
of diffuse emission measu red by RXTE is rather broad 
l|Valinia fc Marshall Il998h so the use of this spectrum is 
appropriate 

A rapid softening of the spectrum below the SPI range 
is indic ated by a comparis on with XMM and Chandra 
results. lHands et all l)2004 (XMM) find a total 2-10 
keV intensity 9.6 x 10~ n erg cm" 2 s _1 deg~ 2 for / = 
19° - 22°, |6| < 0.5°, 80% of which is diffuse or at least 
not in known sources, corres ponding to 0.2 MeV cm" 2 
s" 1 sr" 1 in the units of Fig 171 Ebisawa et alJ (|200ll 12005^ 
(Chandra) finds a similar value for diffuse emission at 
I — 28°, b — —0.2°. Even allowing for the large concentra- 
tion of the emission to the ridge, this is above the extrap- 
olation of the SPI power-law by at least a factor 3. This is 
consistent with the inclusion of a cutoff power-law at low 



Positronium fraction defined as / = 2/ ( 



|r) where 

r is the ratio of 511 keV line flux to 37 continuum flux 
jBrown fc Leventhal fl987h 



6 The total energy flux in this region in the 3-20 keV band 
is 3.5 x 10 -9 erg cm -2 s _1 (M. Revnivtsev, private communica- 
ti on) and this i s used for normalization of the RXTE spectrum 
of lRevnivtsevI (120031) . 
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energies bv iKinzer et al.l l)l999j) in the OSSE analysis. It 
is also compatible with the steepe r power-law (index > 2 ) 
below 20 keV measured by Ginga ijYamasaki et allll997|b 



galdef ID 47_600203a 



if 10" 



10% 
10 





0.25<l< 9.75 , 350.25<l<359.75 




-9.75<b<-0.25 , 0.25<b< 9.75 


'''v'V | 






I 







1 10 

energy, MeV 



Fig. 7. Comparison of spectrum of diffuse emission (the 
central part of the Galaxy (350° < I < 10°, |6| < 
10°)), with o ther e xperiments. Data from RXTE (cyan) 
llRevnivtsevl l2008h. SPI (blue) (th is work), OSSE (ma- 
genta curves) ijKinzer et all 1 1 999^ (their Table 3 and 
Fig 5, VP 5+16: I = 0, b = 0) , COM PTEL (green 
llStrong et all Il999h and EGRET (red) l|Strong et al. 
l2004b(l . Note that the RXTE data refer to the region 
350° < I < 10°, 2 < \b\ < 10° (see text). 



7.2. Sources: resolved and unresolved 

For 20-100 keV the resolved source/diffuse emission ra- 
tio is about 10 (FigEJ, which is compatible with exist- 
ing estimates of the total source emission in the Galaxy: 
1 Grimm et all l2002h of 2 1 39 erg s" 1 and the to tal dif- 
fuse emission of 10 erg s" 1 llDogiel et alll2002ah . These 
estimates are mainly based on the energy range 2-10 keV 
(largely RXTE-based) so this is only a plausibility check. 



It is of interest to estimate the limitations of the 
present analysis which arise from the SPI instrument it- 
self: the number of sources visible is ultimately limited, 



7 It is important to note the reason for the difference between 
the 'diffuse-dominated' Galactic emission found by XMM and 
Chandra (Hands et al. 2004, Ebisawa et al. 2005) and the 
'source dominated' emission shown in Fig 3. The latter refers 
to a large region of the Galaxy including many strong sources, 
while XMM and Chandra refer to a small field-of-view cho- 
sen to avoid strong sources. Since the emission from the whole 
Galaxy is dominated by the brightest sources, the relative im- 
portance of source and diffuse components depends sensitively 
on the size of region considered. 



at low energies, by source confusion. The SPI angular re- 
sponse has a width of about 2.6°, which means that for 
the inner radian in longitude within 5° of the plane we 
cannot detect more than 60 x 10/2. 6 2 ^100 sources. This 
limits the flux which could ever be explicitly attributed to 
detected sources, as follows. From Tabic 2, using the data 
between 3u and 10er and the first 6 energy ranges, we find 
an average logN(>S)-logS slope of -0.9, i.e. converging to 
low fluxes. There are 57 sources > 3er in 18 - 48 keV in 
the area chosen above, and from this it follows that we 
formally have to go down to 1.6cr to reach the confusion 
limit. Integrating N(S)xS, we find that the total source 
flux increases by about 1.5 going from our adopted 3cr 
(Fig 3) to the confusion limit. This is a very rough esti- 
mate, but illustrates that we are near the limit of what 
could be deduced from the SPI data, and that the total 
flux in detected sources could never be more than 50% 
larger than shown in Fig 3. 

Sources undetectable because of confusion or low flux 
will by definition be interpreted as part of the diffuse emis- 
sion in our analysis. The full analysis must then involve 
detailed modelling of source populations, which will be 
addressed in future work. 



7.3. Diffuse continuum emission 

The broad-band spectrum from 2 keV to 100 GeV, in- 
cluding RXTE, SPI, COMPTEL and EGRET, for 330° < 
I < 30°, \b\ < 5° 8 is compared with predictions us- 
ing the galprop model in Fig |SJ The RXTE spectrum 
is for 350° < I < 10°, 2 < \b\ < 10° so that its in- 
clusion is only indicative, but illustrates the softening of 
the spectr um below 20 keV. The model is the 'optimized' 
one from IStrong et all I 2004bldh . including the cosmic- 
ray generated bremsstrahlung and inverse Compton com- 
ponents. Evidently neither of these processes is suffi- 
cient to account for the observed non-thermal hard X- 
ray spectrum; increasing the inverse Compton by invok- 
ing a steeper electron spectrum leads to an overpredic- 
tion of the EGRET data points. Alternative explana- 
tions of the hard X-ray spectrum have been mentioned 
in the Introduction, but these do not easily produce the 
observed very hard spectrum, so that a population of 
compact sources seems to us most likely, and anoma- 
lous X-ray pulsars (AXPs), with th eir very hard spectra 
l|Kuiner. Hermsen fc Mendez 1 . 12004^ . appear good candi- 
dates. AXPs also have an appropriate young population 
distribution compatible with the SPI result. 

This picture is consistent with a change of origin from 
truly diffuse emission below 20 keV with a steep spec- 
trum, as seen by RXTE, XMM and Chandra, to a source- 
dominated origin with a hard spectrum above 20 keV as 
seen by SPI. 

The present work will form the basis for more detailed 
evaluation of these models. 



region chosen for compatibility with lStrong et all J2004b) 
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Fig. 8. Broadband s pectrum of di f fuse em ision (330° < I < 30°, |6| < 5° ), compared with model based on cosmic- 
ray interactions from lStrong et alJ l|2004bllcl . Data fro m RXTE (cyan), SPI (red lines), IBIS (magenta), COMPTEL 
(green), (references as in Figd and EGRET (red bars) l|Strong et all l2004bj) are shown. Dashed/ green curve: inverse 
Compton, dotted/red: 7r°-decay, dashed-dot/light blue: bremsstrahlung, EB: extragalactic background, dark blue: total 
model. Notes: (a) the RXTE data refer to the region 350° < I < 10°, 2 < \b\ < 10° , so this may explain in part the 
high intensity rel ative to the SPI results: (b) the EGRET data >10 GeV have been corrected for the new sensitivit y 
estimate given bv iThompson. Bertsch fc O'Neal (2005), which results in higher values than in lStrong et al ] ll2004bh : 
a 50% flux error has been adopted to plot this data. 
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8. Conclusions 

The diffuse continuum emission from the Galactic plane 
has been detected at high significance with SPI. The 
positronium and non-thermal components are clearly sep- 
arated, with the positronium dominating above 300 keV in 
the inner Galaxy. The contribution from detectable point- 
sources has for the first time been explicitly accounted for, 
which was not possible with earlier experiments. The com- 
bined emission from detectable sources varies from a factor 
10 times the diffuse emission at low energies to negligible 
at high energies. 

Spectral analysis, using the SPI energy response, re- 
veals a power-law component with index about -1.7, con- 
sistent with OSSE and COMPTEL results. The positro- 
nium flux drops by factor at least 7 from |^| < 10° to 
|Z| > 10° in accord with other studies. The power-law 
component by contrast drops by only a factor 2, proving 
a quite different longitude distribution and spatial origin. 

We suggest a change of origin from truly diffuse emis- 
sion below 20 keV with a steep spectrum, to a source- 
dominated origin with a hard spectrum above 20 keV. 

From the broadband spectrum of the diffuse emision 
from keV to TeV energies we see that the emission in the 
SPI range cannot be easily explained as inverse Compton 



or bremsstrahlung from cosmic-ray electrons. Either an- 
other mechanism is at work or there is a new population 
of compact sources with characteristic high-energy tails. 

In view of the apparent difficulty to produce the 
very hard spectrum with a diffuse mechanism, at present 
the source population seems the most plausible explana- 
tion. Possible candidates are the anomalous X-ray pulsars 
(AXPs) with the characteristic hardening of their spec- 
tra around 100 keV. Detailed evaluation of the contri- 
bution from unresolved point-sources will require popu- 
lation synthesis studies, and this will be the subject of 
future work. Meanwhile we are investigating further the 
possibility o f diffuse mechanisms like that proposed by 
iDogiel et all ll2002bl) . 



Appendix A: Model Fitting Details 

The objective of the model fitting analysis is now to ex- 
tract information about 0i in the form of posterior prob- 
ability distributions, and their moments (mean, standard 
deviation etc) and any other functions of interest (e.g. the 
total image). 



10 



Strong A.W. et al.: Gamma-ray continuum emission from the inner Galaxy 



Fitting and error estimation 

The likelihood function is: 



L(D\9)=He- d "d n k "/ 



no- 



where rik are the measured data (denoted collectively by 
D). 

and the posterior probability P(9\D) is expressed in 
terms of the likelihood and the prior probability Pr(0) 
using Bayes theorem: 

P(9\D) = L(D\9)Pr(9)/P(D) 

The posterior for one parameter is obtained by 
marginalizing over the other parameters: 

P(9 i \D)= f P(9' l \D)d N 0' 
and its mean value is 



<9 t \D>= J 9 l P{9 l \D)d9 l = J 9 l P(9\D)d N 9 
with standard deviation 

A9i\D = sqrt ((9- < 9 Z \D >) 2 P(9 l \D)d9 l 



sqrt 



i- < 



\D >fP(9\D)d N 9 



An analytical approximation to the covariance matrix 
and hence the error estimates can be obtained by express- 
ing the log-likelihood function as an expansion about the 
maximum: 



log L(9) = log L(9 ) 



1. 



-£ p X<j( 



Z)(9 q - 9°)H P 



where the Hessian matrix is H pq = d 2 log L/d9 p d9 q . The 
marginalization integr ation can th en be done analytically 
using diagonalization (Sivial Ijl997l) : Ch 3.2 and Appendix 
A. 3) which gives the result for the covariance under the 
posterior: 

a lq =< {Op - ®p)(9q - Oq) >= Hpq 1 

The inverse of H can be written in terms of the matrix of 
eigenvectors X and eigenvalues A = diag(X q ) : 

HX = XA 

X = H^XA 

H- 1 = XA^X- 1 = XA~ 1 X T 

since eigenvectors of a real symmetric matrix are orthog- 
onal, X^ 1 = X T . For a single parameter, 



X 



It is interesting to note that a completely different formal- 
ism i nvolving the d istribution of the maximum-likelihood 
ratio ( Stronel 1 19851) and maximizing the likelihood under 
constraints instead of marginalizing over unwanted param- 
eters, actually leads to the same formula for errors on one 
parameter, once we express the inverse Hessian in terms 
of eigenvalues. This has not been noted before to the au- 
thor's knowledge. One advantage of the Bayes formalism 
is that the covariances can be used to compute the errors 
on linear combinations of components. 

The RMS error on a linear combination of parameters 
is directly related to the covariances as follows 9 

/ = T,w l 9 i 

Af = Wi A9i 
(A/) 2 = ,„;„■, A/>, AO, 



< (A/) 2 >= ZiY.jWiWj < A9iA9.j > 
= Y. l Y. j WiWj<Ji j 



Appendix B: Fit results 

This appendix consists of Tables Bl-4, containing all the 
SPI results plotted in the spectra, to enable convenient 
use of the present work. 
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Table B.l. SPI diffuse emission spectrum for various 

source variability timescales. * indicates selected sources 

(Sco X-l, 4U1700-377, OAO 1657-4154) allowed to vary 

on timescale of SPI pointings. 

0.25° < I < 29.75°, 330.25° < I < 359.75°, 

-9.75° < b < -0.25°, 0.25° < b < 9.75° 



Energy range 


variability 


intensity 


error 


ID 


(ke V J 


timescale 
(days) 


cm - sr~ 
s Mcv 






18.0- 28.0 


const 


5.469 


1.711 


2481 


18.0- 28.0 


200 


4.484 


1.580 


2501 


18.0- 28.0 


100 


3.572 


1.452 


2502 


ion oo n 

18.0- 28.0 


50 


2.870 


1.369 


2503 


18.0- 28.0 




4.341 


2.496 


2505 


28.0- 38.0 


const 


3.725 


0.782 


2482 


28.0- 38.0 


200 


3.112 


0.694 


2511 


28.0- 38.0 


100 


2.846 


0.656 


2512 


28.0- 38.0 


50 


2.626 


0.629 


2513 


28.0- 38.0 




3.295 


0.723 


2515 


38.0- 48.0 


const 


1.894 


0.643 


2483 


38.0- 48.0 


200 


1.861 


0.421 


2521 


38.0- 48.0 


100 


1.783 


0.410 


2522 


38.0- 48.0 


50 


1.778 


0.413 


2523 


38.0- 48.0 




1.826 


0.417 


2525 


48.0- 58.0 


const 


1.859 


0.483 


2484 


48.0- 58.0 


200 


1.958 


0.923 


2531 


48.0- 58.0 


100 


1.977 


0.930 


2532 


48.0- 58.0 


50 


1.786 


0.904 


2533 


48.0- 58.0 


* 


1.755 


0.769 


2535 


58.0- 68.0 


const 


1.685 


0.469 


2485 


58.0- 68.0 


200 


1.643 


0.466 


2541 


58.0- 68.0 


100 


1.602 


0.461 


2542 


58.0- 68.0 


50 


1.559 


0.459 


2543 



Table B.2. SPI diffuse emission spectrum for 
330° < I < 30°, \b\ < 10° 



Energy range 
(keV) 


components 


intensity 
cm~ 2 sr _1 
s" 1 MeV" 1 


error 


ID 


18.0- 28.0 


1-9 


5.469 


1.711 


2481 


28.0- 38.0 


1-9 


3.725 


0.782 


2482 


38.0- 48.0 


1-9 


1.894 


0.643 


2483 


48.0- 58.0 


1-9 


1.859 


0.483 


2484 


58.0- 68.0 


1-9 


1.685 


0.469 


2485 


68.0- 78.0 


1-9 


1.077 


0.232 


2486 


78.0- 88.0 


1-9 


0.378 


0.110 


2487 


88.0- 98.0 


1-9 


0.863 


0.220 


2488 


98.0- 138.0 


1-9 


0.444 


0.024 


2489 


138.0- 178.0 


1-9 


0.265 


0.030 


2490 


178.0- 258.0 


1-9 


0.196 


0.020 


2291 


258.0- 338.0 


1-9 


0.134 


0.013 


2248 


338.0- 418.0 


1-9 


0.106 


0.014 


2249 


418.0- 498.0 


1-9 


0.093 


0.016 


2183 


508.0- 514.0 


1-9 


0.465 


0.064 


2207 


518.0-1018.0 


1-9 


0.018 


0.005 


2186 



Table B.3. SPI total source emission spectrum for 
330° < I < 30°, |6| < 10° 



Energy range components intensity error ID 



l Kev ) 




m — 2 sr~ 1 

O 1V1C V 






ion 09. O 
1(5. U- Zo.U 


in 1 QO 


D4.U04 


U.DoU 




28.0- 38.0 


10-132 


18.619 


0.161 


2482 


38.0- 48.0 


10-106 


9.585 


0.095 


2483 


48.0- 58.0 


10-71 


5.004 


0.116 


2484 


58.0- 68.0 


10-53 


2.707 


0.108 


2485 


68.0- 78.0 


10-45 


1.579 


0.036 


2486 


78.0- 88.0 


10-50 


1.198 


0.029 


2487 


88.0- 98.0 


10-39 


0.770 


0.047 


2488 


98.0- 138.0 


10-32 


0.379 


0.012 


2489 


138.0- 178.0 


10-21 


0.114 


0.010 


2490 



Table B.4. SPI diffuse emission spectrum for 

0.25° < I < 9.75°, 350.25° < / < 359.75°, -9.75° < b < 

-0.25°, 0.25° < b < 9.75° . 



Energy range 
(keV) 


components 


intensity 

cm~ 2 sr _1 
s" 1 MeV" 1 


error 


ID 


18.0- 28.0 


1-9 


13.829 


3.434 


2481 


28.0- 38.0 


1-9 


7.015 


1.402 


2482 


38.0- 48.0 


1-9 


4.502 


1.065 


2483 


48.0- 58.0 


1-9 


3.754 


0.858 


2484 


58.0- 68.0 


1-9 


2.335 


0.680 


2485 


68.0- 78.0 


1-9 


2.217 


0.442 


2486 


78.0- 88.0 


1-9 


0.574 


0.173 


2487 


88.0- 98.0 


1-9 


1.599 


0.376 


2488 


98.0- 138.0 


1-9 


0.825 


0.037 


2489 


138.0- 178.0 


1-9 


0.492 


0.043 


2490 


178.0- 258.0 


1-9 


0.359 


0.028 


2291 


258.0- 338.0 


1-9 


0.237 


0.019 


2248 


338.0- 418.0 


1-9 


0.190 


0.019 


2249 


418.0- 498.0 


1-9 


0.190 


0.022 


2183 


508.0- 514.0 


1-9 


1.244 


0.095 


2207 


518.0-1018.0 


1-9 


0.022 


0.008 


2186 
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Table B.5. SPI diffuse emission spectrum for 

9.25° < I < 29.75°, 330.25° < I < 350.75°, -9.75° < b < 

-0.25°, 0.25° < b < 9.75° 



Energy range 


components 


intensity 

cm - ^ sr ~^ 
s -1 MeV" 1 


error 


ID 


lo.U- Zo.U 


1 Q 

i-y 


l.Oul 


i 1 7Q 


949.1 

Z4:Ol 


oo n oon 

ZO.U - OO.U 


1 Q 
i-y 


9 1 A7 

Zr. ±4: ( 


U.OU4: 




oo.U- 4:O.U 


1 Q 
i-y 


U.D1U 


U.4:00 


9/1 Q. 1 } 


4:0. U- Oo.U 


1 Q 


u.yoo 


u.ooy 


9/1 ft/1 


p:o n fic n 
oo.u- uo.u 


1 Q 

i-y 


l.OUO 


U. 4:4:1 


Z4:00 


Oo.U- (O.U 


1 Q 
i-y 


u.ouo 


U. 104: 


94ftfi 

Z4:O0 


/O.U- oo.U 


1 Q 

i-y 


U.ZoD 


U. 1U1 


9/187 

Z4:0 f 


oo n no n 
oo.U- yo.U 


i-y 


U.4y4 


A 1 79 

U.i /o 


O/IOQ 

z4oo 


98.0- 138.0 


1-9 


0.252 


0.027 


2489 


138.0- 178.0 


1-9 


0.152 


0.033 


2490 


178.0- 258.0 


1-9 


0.119 


0.022 


2291 


258.0- 338.0 


1-9 


0.083 


0.015 


2248 


338.0- 418.0 


1-9 


0.064 


0.015 


2249 


418.0- 498.0 


1-9 


0.046 


0.017 


2183 


508.0- 514.0 


1-9 


0.085 


0.072 


2207 


518.0-1018.0 


1-9 


0.016 


0.006 


2186 



